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Abstract 

This paper is concerned with the coherent quantum filtering (CQF) problem, where a quantum observer is cascaded in a 
measurement-free fashion with a linear quantum plant so as to minimize a mean square error of estimating the plant variables 
of interest. Both systems are governed by Markovian Hudson-Parthasarathy quantum stochastic differential equations driven by 
bosonic fields in vacuum state. These quantum dynamics are specified by the Hamiltonians and system-field coupling operators. 
We apply a recently proposed transverse Hamiltonian variational method to the development of first-order necessary conditions 
of optimality for the CQF problem in a larger class of observers. The latter is obtained by perturbing the Hamiltonian and 
system-field coupling operators of a linear coherent quantum observer along linear combinations of unitary Weyl operators, 
whose role here resembles that of the needle variations in the Pontryagin minimum principle. We show that if the observer is 
a stationary point of the performance functional in the class of linear observers, then it is also a stationary point with respect 
to the Weyl variations in the larger class of nonlinear observers. 


I. INTRODUCTION 

The dynamics of a wide class of open quantum systems, interacting with the environment, can be described in the frame¬ 
work of the Hudson-Parthasarathy quantum stochastic calculus [9], [10], [19]. This approach employs quantum stochastic 
differential equations (QSDEs) which are driven by a quantum mechanical analogue of the classical Wiener process [13]. 
The quantum Wiener process models a heat reservoir of external fields and acts on a boson Fock space [21], [19]. The drift 
vector and the dispersion matrix of the QSDEs depend on the system Hamiltonian and the system-field coupling operators. 
These energy operators specify the evolution of the system as a result of its internal dynamics influenced by the interaction 
with the environment. The Hamiltonian and the coupling operators are usually modelled as functions of the system variables. 
In particular, such functions can be polynomials or Weyl quantization integrals [5], [24], which affects the complexity of 
the resulting quantum system. 

An important role in the linear quantum control and filtering theory [3], [12], [22] is played by open quantum harmonic 
oscillators (OQHOs) [4], [6] with quadratic Hamiltonians and linear system-field coupling operators. The linear-quadratic 
dependence of the energy operators on the system variables, in combination with the canonical commutation relations (CCRs) 
between the variables, makes the OQHO dynamics linear (and Gaussian in the case of vacuum fields and Gaussian initial 
states [11], [20]). Despite some similarities to the classical linear SDEs, the coherent (that is, measurement-free) quantum 
counterparts [18], [17] to the classical LQG control and filtering problems [1], [14] for OQHOs are complicated by the 
physical realizability (PR) constraints. The latter are associated with the state-space matrices of the QSDEs for fully quantum 
controllers or filters and are related, in particular, to the CCR preservation. 

We mention one of the existing variational approaches [27], [28] to the coherent quantum LQG (CQLQG) control and 
coherent quantum filtering (CQF) problems which develops optimality conditions using the Frechet differentiation of the 
LQG cost with respect to the state-space matrices. The quantum nature of the underlying problem enters this approach 
only through the PR constraints, with all the other aspects of the method being essentially “classical”. The latter has certain 
advantages, such as practical applicability to the numerical optimization algorithms [25]. However, this approach is limited to 
linear controllers and filters, and the resulting optimality conditions do not provide insights into whether nonlinear quantum 
controllers or filters can outperform the linear ones for linear quantum plants. 

In the present paper, we consider a CQF problem, similar to [17], [28], where a quantum observer is cascaded in a 
measurement-free fashion with a linear quantum plant so as to minimize a mean square error with which the observer 
variables approximate linear combinations of plant variables of interest. Both systems are governed by Markovian Hudson- 
Parthasarathy QSDEs driven by bosonic fields in vacuum state. We employ a recently proposed fully quantum variational 
method of [30] based on using a transverse Hamiltonian, an auxiliary time-varying operator which encodes the propagation of 
perturbations through the unitary system-field evolution. We apply the transverse Hamiltonian approach to the development 
of first-order necessary conditions of optimality for the CQF problem in a larger class of observers. The latter is obtained 
by perturbing the Hamiltonian and system-field coupling operators of a linear coherent quantum observer along linear 
combinations of the unitary Weyl operators [5]. Similar trigonometric polynomials of quantum variables have recently been 
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used in [24] for modelling the uncertainty in system Hamiltonians. In the present paper, the Weyl variations play a different 
role which resembles that of the needle variations in the proof of the Pontryagin minimum principle [23]. We show that if 
the observer is a stationary point of the cost functional in the class of linear observers, then it is also a stationary point with 
respect to the Weyl variations (with the latter leading to nonlinear observers). Therefore, in the mean square optimal CQF 
problem for linear quantum plants, linear coherent quantum observers are locally sufficient at least in the sense of the Weyl 
variations of the energy operators. 

The paper is organised as follows. Section HIl outlines notation used in the paper. Section HII] specifies the class of open 
quantum stochastic plants being considered. Section |IV] formulates the mean square optimal CQF problem. Section |V] 
employs the transverse Hamiltonian approach in order to study the sensitivity of the performance criterion to perturbations 
in the energy operators of observers. Section [VTI applies these results to liner-quadratic perturbations of the observer energy 
operators and obtains the first-order necessary conditions of optimality in the class of linear observers. Section IVlII introduces 
the Weyl variations of the energy operators and establishes the main result of the paper that the stationarity of the quadratic 
cost functional with respect to linear-quadratic perturbations of a linear observer implies the stationarity with respect to the 
Weyl variations. Section IVIIII provides concluding remarks. 

H. PRINCIPAL NOTATION 

In what follows, {A,B] := AB — BA denotes the commutator of linear operators A and B on a common space. As a linear 
superoperator, the commutator with a fixed operator A, is denoted by ad, 4 (-) := [A, •]. This extends to the commutator {n x m)- 
matrix := for a vector of operators Ai,...,A„ and a vector T of operators 

Vectors are organized as columns unless indicated otherwise, and the transpose (•)^ acts on matrices of operators 
as if their entries were scalars. In application to such matrices, (•)^ := ((•)^)^ denotes the transpose of the entry-wise 
operator adjoint {■)^. For complex matrices, is the usual complex conjugate transpose (•)* := ((•))^- The subspaces of 
real symmetric, real antisymmetric and complex Hermitian matrices of order ii are denoted by S„, A„ and H,, := S„ + iA„, 
respectively, where i := i/—1 is the imaginary unit. The real and imaginary parts of a complex matrix are denoted by Re( ) 
and Im(-). These extend to matrices M with operator-valued entries as ReM = j(M +M*) and ImM = j-.(M — M**) which 
consist of self-adjoint operators. Also, S(M) := j{M + M'^) denotes the symmetrizer of square matrices. Positive (semi-) 
definiteness of matrices and the corresponding partial ordering are denoted by ()^) Also, §+ and EI+ denote the sets of 
positive semi-definite real symmetric and complex Hermitian matrices of order n, respectively. The tensor product of spaces 
or operators (in particular, the Kronecker product of matrices) is denoted by 0. The tensor product A 0 B of operators A and 
B acting on different spaces will sometimes be abbreviated as AB. The identity matrix of order n is denoted by /„, while 
the identity operator on a space H is denoted by J^//. The Frobenius inner product of real or complex matrices is denoted 
by {M,N) := Tr(M*A). Also, ||v||a: := V\AKv denotes the Euclidean (semi-)norm of a real vector v associated with a real 
positive (semi-)definite symmetric matrix K. The expectation := Tr(pi*) of a quantum variable ^ over a density operator 
p extends entrywise to vectors and matrices of such variables. 

HI. QUANTUM PLANTS BEING CONSIDERED 

We consider a quantum plant which is modelled as a quantum stochastic system interacting with m external boson fields. 
The plant has n dynamic variables Ai (f),... ,A„(f) which evolve in time f ^ 0. The plant variables are self-adjoint operators 
on a composite plant-field Hilbert space where is the initial complex separable Hilbert space of the plant which 

provides a domain for Ai(0),...,A„(0), and is a boson Eock space [19] for the action of quantum Wiener processes 
VTi (f),..., Wm(f). The latter are self-adjoint operators which model the external boson fields. The energetics of the plant- 
field interaction is specified by the plant Hamiltonian H{t) and the plant-field coupling operators Li (f),... which are 

self-adjoint operators, representable as time-invariant functions (for example, polynomials with constant coefficients or Weyl 
quantization integrals [5]) of the plant variables Ai(f),... ,A„(f). Therefore, both//(O) and Li(0),... ,Lm(0) act on the initial 
space Omitting the time arguments, we assemble the plant and field variables, and the coupling operators into vectors; 


Ai' 


'Wi' 


'Li 


, W:= 


, L;= 


_A„_ 


Wm_ 


J-'m_ 


Eor what follows, the plant is assumed to be an OQHO [4]. More precisely, the plant variables satisfy the Weyl CCRs 
qi{u+v) X _ ^iu &v^tu Xgiv X fQj. y g where e“' ^ is the unitary Weyl operator [5] on ,^0,^ which is associated 
with (and inherits time dependence from) the plant variables. The Heisenberg infinitesimal form of the Weyl CCRs is 
described by 

[X,X'^]=2i@ (2) 
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where the CCR matrix 0 G A„ represents 0 (g) and remains unchanged. The Hamiltonian of the OQHO is a quadratic 

polynomial of the plant variables, and the plant-field coupling operators in ([TJ are linear functions of the variables: 


1 T 

H ^ -X^RX, 

2 


L = NX, 


(3) 


where R G is the energy matrix, and N G is the coupling matrix. The plant and the external fields form an isolated 
quantum system whose evolution is described by a unitary operator U{t) on driven by the fields and their interaction 

with the plant according to the QSDE [10], [19] 

- (^i{Hodt+LldW) + ^Llaiodt^U 
1 


dU = -{ 


-u(^i{Hdt+L^dW)^ 


-L^Q.Ldt . 


(4) 


with initial condition Uq := The subscript (•)o indicates the initial values of time-varying operators (or vectors and 

matrices thereof), so that Hq :=H{Q), Lq :=L{Q) and Uq := U{Q), and the time arguments will often be omitted for brevity. 
Due to the continuous tensor product structure of the Fock space [21], the future-pointing increments dW commute with 
adapted processes (including U) taken at the same (or an earlier) moment of time. The matrix £2 S IHI+ in 

(|4]i is the Ito matrix of the quantum Wiener process W: 


and the matrix J G 


dWdW * = £idf, a :=!„, + iJ, 
specifies the cross-commutations between the entries Wi 

[dW,dW'^]=2iJdt, 7:=J04/2, J: 


0 

-1 


of W: 

1 
0 


(5) 


(6) 


where the dimension m is assumed to be even. The QSDE ([Hi corresponds to an important particular case of open quantum 
dynamics when the scattering matrix is the identity matrix, and there is no photon exchange between the fields, thus 
eliminating the gauge processes [19] from consideration. A plant operator Ob on the initial space (which can be identified 
with its extension to the plant-field space M’ ®^') evolves to an operator (7(f) on ^ at time f ^ 0 according 

to the flow 


a(t) :=Jr((Jo) = 


(7) 


When it is applied to vectors and matrices of operators, the flow Jf acts entrywise. In view of the identity {gq® = t/<7, 

which follows from (j?]) and the unitarity of U, the second equality in (HI employs the representation of the Hamiltonian and 
the coupling operators in terms of the flow: 

//(f)=j,(//o), L(f)=j,(Lo). 

Note that the flow J; depends on the energy operators Hq and Lq (or the energy and coupling matrices R and A in Q in 
the case of OQHOs). More precisely, any perturbation of Hq and Lq (as functions of the fixed set of system variables Xq) 
modifles the flow. Now, the quantum adapted process a in (|7]i satisfles the following Hudson-Parthasarathy QSDE [10], 
[19]: 

da = ^{a)dt-i[a,L^\dW, <:^{a)-.= i[H,a] + Si{a). (8) 

Here, ^ is the Gorini-Kossakowski-Sudarshan-Lindblad (GKSL) decoherence superoperator [7], [15] which acts on a as 

&{a) := ]^{L}Q.[a,L] + [L^,a]Q.L) 


=-[a, L^]Q.L-- [L^Q.L, a] 

= - [(7,a'^]a'^£ 2AA - i cr]. 


(9) 


The last two equalities in (|9]l are convenient for the entrywise evaluation of ^ at vectors of operators. The superoperator 
in ® is referred to as the GKSL generator. In application to the vector X of plant variables, the flow (|2l) acts entrywise as 

X{t) :=i,iXQ) = U{t)\XQ®J^)U{t). 

In view of © and Q, the corresponding QSDE ([8]) takes the form 

dX = <^{X)dt - i[X,LJ]dW = AXdt -f BdW, (10) 

with the «-dimensional drift vector ^{X) =AX and the dispersion {n x fM)-matrix —i\X,lJ] = — =B, where the 
matrices A G R"^" and B G R"^"' are given by 


A :=2&(R + N'^JN), 


B := 20A' 


( 11 ) 
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The interaction of the input field with the plant produces an m-dimensional output field 

'Yiity 

Y{t):= : =t/(f)i-(^^(g)lT(f))t/(f), 


( 12 ) 


Mt)_ 

where the plant-field unitary evolution is applied to the current input field variables (which is closely related to the innovation 
role of the quantum Wiener process W in the QSDEs). The output field satisfies the QSDE 


dT = IJLdt YdW ^CXdt + m, (13) 

where J is the matrix from (|6]l, L is the vector of plant-field coupling operators in ([Til and (|3ll, and the matrix C € is 

given by 

C:=2JN. (14) 


IV. COHERENT QUANTUM EILTERING PROBLEM 

Consider a measurement-free cascade connection of the quantum plant, described in Section|III] with another open quantum 
system. The latter plays the role of a coherent quantum observer and is driven by HE (which is part of the plant output Y 
in (fT2li ') and a quantum Wiener process co of even dimension /i on a boson Lock space 5^; see Eig. [T] The matrix If £ 



plant 




HE 


observer 





IV 


(0 


Fig. 1. The cascade connection of a quantum observer with a quantum plant, mediated by the field flF and affected by the environment through the 
quantum Wiener processes W and co. Also shown is the observer output 7]. 


is formed from conjugate pairs of rows of a permutation matrix of order m, with p ^m. The observer has its own initial 
Hilbert space f), dynamic variables ,..., i*v with a CCR matrix r? G Ay, and a {p + /i)-dimensional output field 77 : 


(15) 


We denote the observer Hamiltonian by E, while the vectors of operators of coupling of the observer with the selected plant 
output YIY and the quantum Wiener process (O are denoted by 



■^r 


■®r 


' li ' 



, (o:= 


, V ■= 









0:= 


'or 


■‘pr 


ji 




.'Yp. 


(16) 


respectively. The Hamiltonian F and the coupling operators <I>i,..., Op and ‘Pi,..., 'Pp are functions of the dynamic variables 
^ 1 ,..., of the observer and hence, commute with functions of the plant variables, including the plant Hamiltonian H and 
the plant-field coupling operators in L. The plant and the observer form a composite open quantum stochastic system, whose 
vector ^ of dynamic variables satisfies the CCRs 


[jr, =2i0, 


0 : = 


© 0 
0 t? 


^:= 


X 


(17) 


and is driven by a combined quantum Wiener process W with the Ito table 

Q. 0 

0 13 


fl .= 


w-.= 


(18) 


= flAt, 

Here, U is the Ito matrix of the quantum Wiener process co of the observer which is defined similarly to H in (|5ll and (|6]i: 

dcudo)^ = Udf, U :=/p-|-;j0/p/2- (19) 

The Hamiltonian H of the plant-observer system and the vector L of operators of coupling with W can be computed by 
using the quantum feedback network formalism [8] as 

rL4-nT<i>' 


H = r-i-o‘n7L, 


‘P 


( 20 ) 
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While the plant dynamics, governed by (fTOl i and (fOl l. remains unaffected by the observer, the dynamic variables of the latter 
in (flSl) are governed by the QSDE 

d^=G{^)dt-i[^,L^]dW 

= (i[r,^]+A(^))dr-/[^,<i>T]ndy-/[^,'ET]d(u. (2i) 

Here, 

G(C):=i[H,C]+D(C) (22) 

denotes the GKSL generator for the plant-observer system, and 

D(C) = -[C,LT]f2L-i[LTf2L,C] (23) 

is the corresponding decoherence superoperator, similar to (|9|l. In d^lT i. use is also made of the partial GKSL decoherence 
superoperator A which acts on the observer variables as 

A(^) =- - 

- ^ ^ ] (24) 

in view of (fT7l i- (l20l i. with finil^ € H+ being the quantum Ito matrix of YIW. Now, similarly to [17], [28], we formulate a 
CQF problem as the minimization of the steady-state mean square discrepancy 

^:=limEZ(f)—^ min, Z:=E^E, E:=FX-G^ (25) 

between q linear combinations of the plant variables of interest and observer variables as specihed by given matrices F G 
and G G (with E being interpreted as the estimation error). Here, the quantum expectation E(-) is taken over the tensor 
product p := 07 0 V of the initial quantum state 07 of the plant-observer system and the vacuum state v in the composite 
boson Lock space ^0^ for the external helds. Also, it is assumed that the plant-observer system is ergodic, whereby the 
limit in (i25l l reduces to averaging over the invariant state of the system, provided the plant and observer variables satisfy 
an appropriate version of the uniform mean square integrability condition. The criterion process Z in (i25T l is a quadratic 
function of the vector JZ" from (Eli: 

Z := ^ := [F -G]. (26) 

The minimization in (i25l l is carried out over the observer Hamiltonian T and the vector <I> of the observer-plant coupling 
operators in (0, while 'T and all the dimensions are hxed. This problem extends [17], [28] in that we do not restrict 
attention to linear observers even though the plant is an OQHO. 


V. INFINITESIMAL PERTURBATION ANALYSIS 

In order to develop hrst-order necessary conditions of optimality for the CQF problem (i25l l in an extended class of 
observers, we will apply the transverse Hamiltonian variational method of [30] to the inhnitesimal perturbation analysis of 
the performance criterion To this end, suppose To and <I>o depend smoothly (for example, linearly) on a small scalar 
parameter e and are perturbed in the directions 

Ko := Eg, Mg := Og, (27) 


consisting of self-adjoint operators on the observer initial space Jo, representable as functions of the observer variables, with 
(•)' := <3£(-)|e=o- The corresponding perturbations of the plant-observer Hamiltonian and the coupling operators in (i20l) are 


Hg — Kg -f Mg nyLg, 


K = 


n^Mg' 

0 


(28) 


The propagation of these perturbations of the energy operators through the unitary evolution U of the plant-observer-field 
system on the space 0S^0 ^ 0^ is encoded by the transverse Hamiltonian [30]. The latter is a time-varying self-adjoint 
operator dehned by Q := /U^U', so that U(f)' = —i\d{t)Q{t) for all f ^ 0, with zero initial condition Qq = 0. In view of (fTSIl . 
dZTl i and ( i28l l. the general QSDE obtained in [30, Theorem 1] for the transverse Hamiltonian, takes the form: 


dQ=[K + hfYUL- Im (y fl )df -f [M^H 0\dW 

= {K-lm{{2L + n^^yQn^M))dt+M'^ndW. 


(29) 


Here, K := jr(Kg) and M := j((Mg) are the evolved versions of the initial perturbations from ( |27] | under the unperturbed 
flow jt{Q ■= (Si C (S {t) of the plant-observer-field system (which is applied here to observer operators 

^ on Jo). In ( |29] |, use is also made of the relation Im(L^£2n^M) = —M^YIJL which follows from (|5]), the commutativity 
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[L,M^] = 0 and the antisymmetry of J. Since the criterion process Z in (l25l l does not depend explicitly on the energy 
operators, the corresponding formal Gateaux derivative of the cost functional can be computed by using [30, Theorem 
2, Section VII] as 

:=^ntn E0(f), 0:=/[2,Z]. (30) 

Here, 0 is the derivative process [30] associated with Z. Its expectation satisfies the integro-differential equation 

(E0)-=/E[e,G(Z)]+E;t(Z), (31) 

where G is the unperturbed plant-observer GKSL generator given by ( i20l i. (i22l i. ( 12^ . and % is an auxiliary linear superoperator 
acting on plant-observer system operators a as 

x{cy) :=/[/:-Im((2L-l-n^O)'^nnTM), a\ 

-2Re([(7, (L + n'^Oj'^jnn'^M). (32) 

Note that z(o’) depends linearly on the perturbations K and M. Now, in addition to the plant being an OQHO, suppose the 

unperturbed observer is also an OQHO with energy matrix r G Sy and coupling matrices N\ G and N 2 G The 

corresponding observer Hamiltonian T and the coupling operators in (ITbl) are 

<I> = Ni^, 'F=JV2^. (33) 

In this case, in view of (ITSl l and (i24l i. the QSDE (i2Tl i becomes linear: 

= a^dt+ l?idy+ I?2d0), (34) 

where the matrices a G bi G R''^'’ and b 2 G R''^^ are computed as 

a := 2r}(r + JV^njn'^JVi -l-VjlmOV2), (35) 

b2:=2mj. (36) 


Therefore, the plant-observer system is governed by a linear QSDE 


d^ = ^^dt + mW, 


A 

b\C a 


B 

b\ b2 


(37) 


Eor what follows, both matrices A and a are assumed to be Hurwitz, and hence, so is This implies that the plant-observer 
system is ergodic and has a unique invariant state which is Gaussian [20] with zero mean and quantum covariance matrix 
^3^ + 10 G Here, G S^+y is the controllability Gramian [1], [14] of the pair ,3§) which is a unique solution of 

the ALE 

sy^ + 3^sy'^ + +S+f = 0. (38) 


Lemma 1: Suppose the plant and the unperturbed observer are OQHOs described by Q, ([T0l i- ([l4l) and (l33T l- (l36l l. with 
Hurwitz matrices A and a. Also, let the observer Hamiltonian and the observer-plant coupling operators are perturbed 
according to (i27l i. Then the corresponding formal Gateaux derivative of the cost functional in (125]) can be computed as 

J^' = Eh,z(ct), (39) 

Here, £2 G S^+y is the observability Gramian of the pair satisfying the ALE 

+ + ^'^^ = 0, (40) 


where the matrices and are given by (126]) and (]T7] ). Also, E*(-) denotes the quantum expectation over the invariant 
Gaussian state of the plant-observer system (with E* S3 =0 and IS := ReE*(j?r S3"^) found from ([38])). Eurthermore, the 
superoperator x from (]32l ) acts on the operator a in ([39]) as 


x{(y)=i 


fG-Im 


Proof: Since Z in 


-f8Im(jrT^0 
is a quadratic form of S3, then (]30l ) 


2V' 

■ AtT ■ 

L^i^n 


nn'M), (7 


nn'^M 


and (1311 ) imply that 




T:=/ hm E[e,E], 

r ^+00 


(41) 

(42) 
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where T G §„+v is the limit mean value of the derivative process /[2,E] associated with E := Re(^^^) = — i0. 

The GKSL generator G of the unperturbed linear plant-observer system acts on E as G(E) = jz/E + Ejz/^ + and 

hence, similarly to [30, Example 3, Section Vll], the matrix T satisfies the ALE 

£/r+r£/'^+E,x(^) = o, (43) 


where the superoperator x in (l32l i is applied to E entrywise. By combining the integral representations of the solutions of 
the ALEs (l40l i. (|4^ and using duality, it follows from (l42l i that 

= (^ y"^“’e*^VT're'-^df, E*;t(E)^ 

= (^,E*x(E)) = E*x(jrT^jr), 


which establishes ( |39] | since ;|f(E) = xi'^ Indeed, a enters xi'^) in (l32l i only through the commutators [d,-] and 
hence, ((7+ 0) = Xi'^) for any 0 G C. The representation (HTI) is obtained by substituting L from Q, <I>, 'E from (l33]) and 
(7 from (l39l l into (l32l i. and using the relation ^[cr, = 20^3^. In particular, the second line of (HTT) is established by 


-2Re([c7, (L + n'^4>)^]f2n'^M) 


= -2Re 






A' 


nn'M 


= -2Re((-4/0^jr)'^ 


A* 


= 81m 


A' 

lA^^H 

nn'^M 




Although the unperturbed observer in Lemma[i|is an OQHO, the perturbations dZTl i are not assumed to be linear-quadratic. 
Therefore, the lemma provides a perturbative tool to develop conditions for such an observer to be a stationary point of the 
CQE problem (l25l l in a wider class of observers. 


VI. OPTIMALITY AMONG LINEAR OBSERVERS 

We will now apply Lemma [T] to the first-order necessary conditions of optimality among linear observers. To this end, 
associated with the Gramians and ^ from ([38]) and ( i40l ) is the Hankelian 


:= 


(fii 

<^21 


^12 

<^22 




(44) 


which, together with 3^ and cS, is split into appropriately dimensioned blocks according to the partitioning of X in 
(O, with (•),■, the jth block-row and (•),^. the kth block-column. 

Theorem 1: Under the assumptions of Lemma [T| the linear observer, described in the lemma, is a stationary point of the 
CQE problem (i25l l in the class of linear observers if and only if 


3^22 G Av, (45) 

n(c4\ +B^^n + blB22)-& = IlJbj322. (46) 

Proof: Let the matrices r and Ai in (1331) be smooth functions of the parameter e, and hence, the corresponding 

perturbations of the linear observer in (IZTl i take the form 

M = a;^, (47) 

where the matrices r' G Sy and A[ G can be arbitrary. Substitution of (i47l i into (HTt leads to 


X{(y) 


2N^ 

NjYl 


nn^A! 




-81m( ^'^^0 


■ A^ ' 
Aj^n 


nn’^A! 


( 5 ). 


(48) 
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where a is the quantum variable from (1^ . The identity = 2^^{R2&Ri —Ri&R 2 )^, which holds 

for any R\,R 2 S S„+v due to the CCRs (fTTl i. implies that 


i[^V^,(T]=2jrT(^0|o 0 

= 2jrT(^[0 ^.2t?r']- 


0 0 

0 r' 

0 

r't>^2. 


0^j JT 
SC. 


The averaging of the latter representation over the invariant quantum state leads to 


^E4^V^,a]=2([0 ^.2t?r'] 



= -4(tJ^ 2 .^. 2 , r') = -4(S(t?<r22), r') , 


where the equality ^ 2 . ■^•2 = <^22 follows from (l44l i. and the symmetry of r' is used. Therefore, in view of ( |39] |. the 
corresponding formal Frechet derivative drSC = —4S(tJ(^22) of the cost functional vanishes if and only if the matrix <^22 
satisfies By a similar reasoning, in view of ( |48] |. the Gateaux derivative of SC along M in Wt\ takes the form 


8ImE*Lr‘^0 


— ReE* 


JT* 


■ AtT ■ 
NjU 

'2N'^' 

NjYl 


Q.n^N\ 


!«) 


(7 


=8Im<^^0 


4Re 


■ AtT ■ 

r2A?T 

nJu 


SlYfN[,3^,2 + 


nn^A?; ,/(#'^0-0^)) 


=8{|ny[A? u^Ni]0S,2-u[n u^Ni]e^,2^ 
+ ujN{sJi t? - &S 12 ) - njn^NiS{iSS22), n[ ^ 

=4 (n((c4\ +B'^^n + bJJS22)^ - JbJS22), N[), 


(49) 


where use is also made of (fTTT) . (fT4li . (l36l) and (l44b . Here, S(tJ<f 22 ) = 0 under the condition (|45]) . in which case, the formal 
Frechet derivative = 4-Il{{CS'2i+B^S2\2+bJS222)'& — JbJ<S 22 ) vanishes if and only if (|46] | holds. Therefore, (|45]) and 
(l46l l are indeed equivalent to the stationarity in the class of linear observers. ■ 


VII. WEYL VARIATIONS OF OBSERVERS 

Now, consider the following bounded perturbations of the observer Hamiltonian and the observer-plant coupling operators 
in (| 22 l); 

K :=Re(aW„) = ^{a\Nu + a\N^u) = |a|cos(M'^^ -f 7 ), (50) 

M:=Re(j3W„) = l(j3W„+^W_„) = {\l5k\co&{u'^X + (51) 

with parameters a € C, )3 := £ C^, 7 ;= arga, 9k ;= argjS^ and u £ M''. Here, W„ ;= = Wi„ is the unitary 

Weyl operator associated with the observer variables. We will refer to the perturbations (l50l) and (fSTI) as the Weyl variations 
(of the observer energy operators). Similar trigonometric polynomials of quantum variables have recently been used in [24] 
to model uncertainties in system Hamiltonians. However, in what follows, the Weyl variations play a different role which 
resembles that of the needle variations in the proof of the Pontryagin minimum principle [23]. 

Theorem 2: Under the assumptions of Lemma [T] the linear observer satisfies the first-order necessary conditions of 
optimality in the CQE problem (i25T l among linear observers if and only if it is a stationary point with respect to the Weyl 
variations in dSOl i and (iSTl i for any a £ C, j3 £ C'’, u £ 

Proof: The stationarity of the cost functional SC with respect to arbitrary Weyl variations of the observer implies the 
stationarity with respect to the linear-quadratic perturbations (i47] | because and = —Re(^^^). In 

fact, arbitrary polynomial perturbations can be reproduced by the Weyl variations which form a larger class of perturbations. 
Therefore, it remains to prove that (l45T l and (i46T l imply the stationarity with respect to arbitrary Weyl variations. By Lemma[T] 
the linear response of the cost functional SC to the Weyl variation K of the observer Hamiltonian in (fSOl l takes the form 

/E 4 Re(aW„),( 7 ] =-Im(aE4W„,(j]). 


(52) 
























Let S[, denote a linear superoperator, which is parameterized by m G and carries out a unitary similarity transformation 


of quantum variables as S„(C) (i^) = (C)> where use is made of a well-known identity for 

operator exponentials [16], [31]. Its application to from (ITtI i leads to S„(,^) = -f 2 '^u since =0 (whence 


S„(X) =X) and = —2/i?m (whereby S„(^) = ^ +2§u). Therefore, in view of (O, S„(a): 

= a — + ^ 22 i^u), and hence, 

[W„,(7] = (S„((7)-(7)W„ 

= -Au^^J22iX + ^22(-? + t?M))W„ 

= 4iu^ ■&i^2id„ + 


^=0' 


(53) 


Here, e‘'^^^\Nu = in view of the commutativity [e'**^^, W„] = 0 between the Weyl operators for all w G R", u G R'', 

and hence. 




vf’=0 ’ 


Also, the relation + d-u)\Nu = —idu\Nu follows from 

^w„ = -/a,,Wv|„^oW„ 

= -/a,,(e'"^'^''W„+.,) - iduWu 


(54) 


(55) 


due to the Weyl CCRs WyW„ = e“'^'^'’W„+y which hold for all m,v G R^. By averaging (f53]) over the invariant Gaussian 


quantum state, whose quasi-characteristic function [2] is E*e 

E*[W„, a] = 


= e^2" 
1 1 


for any A := 


G R"+'', it follows that 


= -4!Vr>^2.^ 


H’=0 

-lhf^22 


(56) 


=-4/MTs(jJ^22)Me 2 "“11^22, 

where use is made of the Hankelian from (i44l i. By substituting (fS^ into (i52l i. the Gateaux derivative of S’ along K in 

takes the form iE*[/r,a] = 4M^S(t?(o22)Me 2lf"ll^22Rea. Therefore, the fulhllment of (145]) makes this derivative vanish for 
any a G C and m G R''. Similar arguments (see also ( i49] l) lead to the following Gateaux derivative of S’ along M in dSB: 


8ImE*(jrT^0 


A/'n 


f2nTRe(]3W„ 


-ReE* 




2N^ 

fA^’n 


nnTRe(j3w„), (T 


=8Im^^0 


A' 


+ Re 


Aj^n 
2A'^ 


f2n‘,E*(jrRe(j3‘W„)) 


Aj^n 


nn\E4(j, jrRe(j3‘w„)] 


= _ 8 e 2l!"llk2lm(^^0 


A* 

A^n 


nn 




+ 4e 2li"llK2Re 


2A^ 

A?’n 


1(0^2 - ^2i^)uimi5'' 


=4lml5'^U{JbJS’22 - {CS’^i +B'^Si2 + bJS22)^H 


(57) 
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Here, Eh,[(T,Wj,] = 0 follows from (l56l l under the condition (|45]) . Also, in view of (l5^ - (l55l l. use is made of the relations 


E*(jrw„) 


E4(j, jrw„] 


—idiEtS 




1 

■ 0 ■ 

1 w=0 

tXu 


E*W„ 






-idxE4a,e‘ 




1 

■ 0 ■ 

lvv’=0 

■du 


EJ(7,W„] 


>^=0 


45;^ (A'^S(0<f )Ae^2ll^ll"^) 
4e^3l!“ll"^22(0#i2-4\i^)"- 


The fulfillment of ( 1451 ) and ( l46b makes the Gateaux derivative ( fSTl i of ^ along the Weyl variation M in ( fSTT i also vanish for 
any j3 € G’ and m G R''. Therefore, (|45] | and (l46l) indeed imply the stationarity with respect to the arbitrary Weyl variations 
® and (EB- ■ 

Note that the proof of Theorem |2] (which employs the Gaussian averaging of quasi-polynomials of system variables, such 
as is closely related to the integro-differential identities [29] for expectations of Weyl quantization integrals over 

Gaussian states. 


VIII. CONCLUSION 

We have applied the transverse Hamiltonian variational method [30] to the mean square optimal CQF problem for a 
linear observer cascaded with a linear quantum plant. It has been shown that if such an observer is a stationary point of 
the problem among linear observers, then it also satisfies the first-order necessary conditions of optimality with respect to a 
wider class of Weyl variations of the energy operators of the observer. In this sense, linear observers are locally sufficient 
for linear quantum plants as far as the mean square performance criteria are concerned. Similar ideas can be developed 
for the CQLQG control problem [18] and an optimal control theory for classical port-Hamiltonian systems [26] and their 
stochastic versions. 
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